Fluctuation dissipation relations in stationary states of interacting Brownian particles under shear 
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The fluctuation dissipation theorem (FDT) is studied close to the glass transition in colloidal suspensions 
under steady shear. Shear breaks detailed balance in the many-particle Smoluchowski equation, and gives 
response functions in the stationary state which are smaller at long times than estimated from the equilibrium 
FDT. During the final shear-driven decay, an asymptotically constant relation connects response and fluctuations, 
restoring the form of the FDT with, however, a ratio different from the equilibrium one. 

PACS numbers: 82.70.Dd, 64.70.P-. 05.70.Ln, 83.60.Df 



In thermal equilibrium, the response of a system to a small 
external perturbation follows directly from thermal fluctu- 
ations of the unperturbed system. This connection is the 
essence of the fluctuation dissipation theorem (FDT) which 
lies at the heart of linear response theory. In non-equilibrium 
systems, much work is devoted to understanding the general 
relation between fluctuation (correlator C{t)) and response 
(susceptibility x(0) functions. It has been characterized by 
the fluctuation dissipation ratio (FDR) X{t) defined via 

It is unity close to equilibrium, X'^'^'>{t) = 1, but deviates in 
non-equilibrium because the external perturbations act against 
non-vanishing currents (see Eq. (|2]i below); FDRs quantify 
the currents and signal non-equilibrium 1 1]. 

Colloidal dispersions exhibit slow cooperative dynamics at 
high concentrations and form glasses. These metastable soft 
solids can easily be driven into stationary states far from equi- 
librium by shearing with already modest flow rates. Spin 
glasses driven by nonconservative forces werepredicted to ex- 
hibit nontrivial FDRs in mean field models 12(1. Such behav- 
ior was observed in detailed computer simulations of sheared 
super-cooled liquids by Berthier and Barrat . During the 
shear induced relaxation, the FDR for particle motion perpen- 
dicular to the shear plane is different from unity, but constant 
in time. This ratio was also found to be independent of ob- 
servable, which led to the notion of an effective temperature 
X = T/Tcs describing the non-equilibrium state. Further 
simulations with shear also saw Toff > T S H I3l re- 
cently Toff was connected to barrier crossing rates [811. 

On the theoretical side, much effort has been made on dif- 
ferent spin-models, close to criticality. Universal FDRs were 
found under coarsening and under shear ifioll . where at the 
critical temperature, a universal value of X = ^ was found. 
Yet, the situation for structural glasses has not been clarified. 

In this letter, we investigate FDT for colloidal suspensions 
close to a glass transition under steady shear starting from the 

-particle Smoluchowski equation. While time dependent 
correlation functions are calculated in the integration through 
transients (ITT) approach 111, T2l, which is based on mode 
coupling theory, the connection to the corresponding suscep- 
tibilities will be derived for the first time here. We show that 



equilibrium FDT is violated, but can be restored in a well de- 
fined sense with a renormalized FDR at long times; however, 
the ratio depends on variable, contradicting the notion of an 
effective temperature. Moreover, we establish a connection 
to the concept of a yield stress, which gives a scenario quite 
different from mean field spin glass (3l. 

N spherical Brownian particles of diameter d, with bare dif- 
fusivity Dq, and interacting via internal forces = —diU, 
i = 1, . . . , iV, are dispersed in a solvent with a steady and ho- 
mogeneous velocity profile v(r) = n r, with shear rate tensor 
K = 7xy. Neglecting hydrodynamic interactions, the distri- 
bution of particle positions evolves according to the Smolu- 
chowski equation llSll 

dt^it)^n-^{t), n = Y,9^-[d,-F,-K-r,], (1) 

i 

where is the Smoluchowski operator and we have intro- 
duced dimensionless units for length, energy and time, d = 
ksT = Do = 1. The Smoluchowski operator for the system 
without shear (n = 0) and the flow-part will be denoted 
and Sn = U — ^le- We distinguish two time-independent dis- 
tributions, Hg^'e = without shear and il'i's ~ for the sta- 
tionary system. Averages are (...) and {. . .)^^\ respectively. 

Stationary correlation functions are Cab (t) = {Sa*e^^''*'5b)^^\ 
where fi^ is the adjoint operator obtained by partial integra- 
tions iflHITill : a fluctuation equals Sa = a — {a)^^\ Note that 
shear in Eq. (|T) leads to a non-Hermitian eigenvalue problem 
1 1511 . The susceptibiUty Xah{t) describes the linear response 
of the stationary expectation value of b to an external pertur- 
bation he{t) shifting the internal energy U toU — a* he{t). 



it) - (6)^^) = / dt'x.b{t - t')K{t') + 0[hl) 



One finds Xab(0 = (E, Itt ' ^^e^ 'b)^^^ M- In non- 
equilibrium, where detailed balance is broken and a nonzero 
stationary probability current jf = [—di + F, + k ■ r^J^'s ex- 
ists, the equilibrium FDT is extended (with j\ the adjoint of 
the current operator defined by jf = ji^ s). 



AXab(t) -XabW + CabW 



-(£i.'£e-\)i^^ (2) 



2 



and a deviation of the fluctuation dissipation ratio (FDR) 

Xab(i) 



-Cab(i) 



(3) 



from unity, the value close to equilibrium, arises. While 



Eq. (O has been known since the work of Agarwal ll4ll . we 
will analyze it for driven metastable (glassy) states and show 
that the additive correction A^ab (t) Ola Il7[ llSll leads to the 
nontrivial multiplicative correction, i.e., a constant FDR at 
long times. For simplicity, we will look at auto-correlations 
(6 = a) of x-independent fluctuations, Sil'^a = 0, where the 
flow-term in the current operator j| in (|2]i vanishes. 

VPs is not known and stationary averages are calculated via 
the ITT approach 111 llll2ll . 



(...)(^) = (...)+7 



ds{lJxy 



where a^y = — ^iVi is a microscopic stress tensor el- 
ement. (Operators act on everything to the right, except for 
when marked differently by bracketing.) ITT simplifies the 
following analysis because averages can now be evaluated in 
equilibrium, while otherwise non-equilibrium forces would be 
required lll9ll . E.g. due to di^e = ^i^e, the expression (|2]i 
vanishes in the equilibrium average. The remaining term is 
split into three pieces containing 17^ 



-7 



2 Jo 

(4) 

We start with the first term in the square brackets (without the 
factor i) which can be integrated over s directly giving 



7 (CT: 



xy 



Otyj 



(5) 



where from now on we consider fluctuations from equilibrium 
5a = a — {a); (the constant (a) cancels in (|4]l ). Intriguingly, 
in Eq. Q the two time correlator enters. 



z^''5a*e^''5a), (6) 



where the rheometer has been shearing for a period before 
the correlation measurement is started. It is one of the central 
quantities in the spin-glass theory of aging While the tran- 
sient correlator cl*^ (t) = Ca(t, 0) = (Ja*e^**(5a} describes 
the dynamics after switch on of the rheometer, the stationary 
correlator Ca(t) = Ca(t, oo) is observed after waiting long 
enough; it measures fluctuations in the stationary state. 

Our approximation for /S.Xa{t) in Eq- © rests on the obser- 
vation that it contains the product of a fluctuation 6a and the 
stationary current. We expect current fluctuations to always 
decay to zero, even in possible non-ergodic situations, and 
thus search for a coupling of Axa(i) to derivatives of Ca(t) as 
they cannot be non-ergodic. Partial integration can be used to 
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FIG. 1: C{t) from the i^^'j"' -model (H and x(i) via Eq. ^ for a 
glassy state (e = 10"'^) and 7 = 10~^" with n — 1...4. Shown are 
integrated correlation, 1 — C{t) and response x'(^) — /J x(J'')dt'. 
Inset shows additionally the normalized transient correlator C'*' for 
comparison and the = i susceptibility for 7 = 10~*. 



the latter term contains the equilibrium derivative fija*. It is 
not conserved and decorrelates quickly as the particles loose 
memory of their initial motion even without shear The lat- 
ter term then is the time derivative of the equilibrium corre- 



lator, Ct'\t) = {Sa*e''^'^Sa). A shear flow switched on at 
t = should make the particles forget their initial motion 
even faster, prompting us to use the approximation e^^ * w 
tp^~nlt^nU ^jjjj projector P = 6a){6a*6a)-'^{Sa* in 



{{nla*)e^''*Sa)\ it is then assured to decay faster than in equi- 
librium. This leads, together with an analogous approxima- 
tion in ((5a*e^^'*(5a), to 



_d_ 

dtu 



c'i\t) 



(7) 



The last term in Q will be identified as short time derivative 
of Ca \ connected with the shear independent decay, while 
■^—Cii{t, tui)\t^=o will turn out to be the long time derivative 

of C'a \ connected with the final shear driven decay. This is 
our main result. It captures the additional dissipation provided 
by the coupling to the stationary probability current in Eq. (|2]). 

In order to proceed, the difference between the stationary 
and the transient correlators needs to be known. We will com- 
ment below on the interesting result for the FDR following 
from the simplest approximation to set them equal. Going be- 
yond this leading approximation can be done via Eq. (|6]l 



Jo 



dty 



(8) 



tu,=0 



show g|;;Ca(i, t^)|t„ 



where we used tw = 00, and factorized the appearing two- 
time average with the projector (Jxy){(TxyCrxy)~^ {c^xy A 
small parameter (J = 7 ds{(TxyC'X-'p{^'^ s)^xy) / {f^xyf^xy) 
C'i\t) - ((17ta*)e^^**5a), where arises which contains as numerator the stationary shear stress 




FIG. 2: Parametric plot of correlation C{t) versus response x'{t) = 
Jg x{'t')dt' for a glassy state (e = 10~^) from the F^'^'-model IH 
together with constant non-trivial FDR (straight lines) at long times. 
The vertical solid line marks the plateau /. Inset shows the FDR 
X{t) as function of strain for the same susceptibilities. 



FIG. 3: Long time FDR as function of shear rate for glassy (e > 0) 
and liquid (e < 0) states in the F^^^-mode\ tz^l , when approaching 
the glass transition for e = ±10^~^'~'^'~*'^ . Inset: lim-y.^o Xq(t 
oo) as function of wavevect or q for incoherent density fluctuations at 
the critical density (e = 0) 12211 . 



measured in 'flow curves' as function of shear rate Ml ill . For 
hard spheres, the instantaneous shear modulus {(JxyCxy) di- 
verges i fisl] giving formally (7 = and that transient and sta- 
tionary correlator agree. In recent simulations of density fluc- 
tuations of soft spheres ||20|] . the difference between the two 
correlators was found to be largest at intermediate times, and 
Ca{t) < C^\t) was observed. Both properties are fulfilled 
by Eq. dill. 

After the discussion of the first term in (|4]i, we turn to the 
correction containing the last two terms in It has van- 
ishing initial value and in a mode coupling approximation in 
ITT for the case of density fluctuations, the two terms also al- 
most cancel each other at long times making their sum a small 
correction. Here, we proceed by ignoring it until a future pre- 
sentation. We hence find 

xa(t) « -a(t) + \ (cl^Ht) - ci^\t)^^\ . (9) 

V Ca [t) J 

In the limit of small shear rates for glassy states, the coiTela- 
tors exhibit two separated relaxation steps ll2l 12111 . During 
the shear independent relaxation onto the plateau of height 



the universal value lim^^o X^{t -> oo) = X(""'^)(7t) = \, 
independent of a. This is in good agreement with the find- 
ings in fsll. The initially additive correction in Eq. (O hence 
turns then into a multiplicative one, which does not depend on 
rescaled time during the complete final relaxation process. 

For a more precise investigation of the FDR, we have to 
consider the difference between the transient and the station- 
ary correlator in Eq. ([8]). We turn to the schematic i^^'^^ -model 



of ITT II22II . which has repeatedly been used to investigate 
the dynamics of quiescent and sheared dispersions 11211 . and 
which provides excellent fits to the flow curves from large 
scale simulations ll23ll . It provides a normalized transient cor- 
relator C^*'' (i), as well as a quiescent one, representing co- 
herent, i.e., collective density fluctuations. The correspond- 
ing stationary correlator C is calculated in a second step via 
Eq. (O. Fig. [T] shows the resulting x together with C for 
a glassy state at different shear rates. For short times, the 
equilibrium FDT is valid, while for long times the suscep- 
tibility is smaller than expected from the equilibrium FDT, 
this deviation is qualitatively similar for the different shear 
rates. For the smallest shear rate, we also plot x calculated by 
(*) 



given by the non-ergodicity parameter /a, we have C^' (t) « Eq. ^ with Ca replaced by Ca, from which the universal 

X(""'^)(7i) = i follows. In the parametric plot (Fig. IS, this 
leads to two perfect lines with slopes — 1 and — i connected by 
a sharp kink at the nonergodicity parameter /. For the other 
(realistic) curves, this kink is smoothed out, but the long time 
part is still well described by a straight line, i.e., the FDR is 
still almost constant during the final relaxation process. We 
predict a non-trivial time-independent FDR Xii{jt) =const. 



ci'^\t), and the equilibrium FDT holds. During the shear- 
induced final relaxation from /a down to zero, i.e., for 7^0, 
and t 00 with t7 = const., the correlator without shear 
stays on the plateau and its derivative is negligible. A non- 
trivial FDR follows. Summarized we find in the glass 



lim Xa(t) 
7^0 




7t < 1, 
^Ci'\t) jt = 0{l). 



It is interesting to note that approximating stationary and tran- 
sient correlator to be equal \12], ci^\t) « Cs,{t), we find 
Xa(i) = — \Ca.{t) for long times. The FDR in this case takes 



if ci*"* (and with Eq. (|8]i also Co) decays exponentially for 
long times, because A^a then decays exponentially with the 
same exponent. The line cuts the FDT line below f for 7 — > 0. 
All these findings are in excellent agreement with the data in 
The FDR itself is of interest also, as function of time (in- 



4 




X=1/2 
incoherent 



10'' 



t [LJ units] 



FIG. 4: Comparison to simulation data for incoherent density fluctu- 
ations in the neutral direction (wavevector q = 7.476^) at tempera- 
ture T = 0.3 (Tc = 0.435) and -y = 10"'^. Circles and squares are 
the data (including units) from Fig. 1 1 in Ref. 0], lines are 1 — Cq 
from Fig. 8 in Ref. (sj], and the response Xq(i) ~ Jq X<i{t')dt' cal- 
culated via Eq. ([9}. The dashed line shows Xq with approximation 



a 



Cq. Inset shows the different correlators, see main text. 



set of Fig.|2]i. A rather sharp transition from 1 to i is observed 
when C^*-' « C is approximated, which ah^eady takes place at 
jt « 10"^, a time when the FDT violation is still invisible 
in Fig. [T] For the realistic curves, this transition happens two 
decades later. The huge difference is strikingly not apparent 
in the parametric plot. 

Fig.[3]shows the long time FDR as a function of shear rate 
for different densities above and below the glass transition, de- 
termined via fits to the parametric plot in the interval [0 : 0.1]. 
In the glass X{t oo) is nonanalytic while it goes to unity in 
the fluid as 7 — > 0, where we verified that the FDT-violation 
starts quadratic in 7 as is to be expected due to symmetries. 
Our analysis also allows to study the interesting question of 
the variable dependence of the FDR, for which we consider 
incoherent, i.e., single particle fluctuations ll22ll which were 
most extensively studied in Jst). The FDR is isotropic in the 
plane perpendicular to the shear direction but not independent 
of wavevector q, contradicting the idea of an effective temper- 
ature as proposed in Q lit] and others (see inset of Fig.|3]l. 

That Eg. ^ is nevertheless not in contradiction to the data 
in Ref. ]l3|] can be seen by direct comparison to their Fig. 11. 
For this, we need the quiescent as well as the transient corre- 



lator as input. (7^*^' has been measured in Ref. ||21[] suggesting 



that it can be approximated by a straight line beginning on the 
plateau of Cq(i). In Fig. |4]we show the resulting suscepti- 



bilities. There is no adjustable parameter, when Cq 
taken, for the other curve, we calculated C^^ (t) from Eq. (O 
using a ~ 0.01 (in LJ units) as fit parameter. The agreement 
is striking. In the inset we show the original Cq from Ref. ^ 
together with our construction of df'^ and the calculated C^^ , 
which appears very reasonable comparing with recent simula- 
tion data on Ca(i, tm) ll20|] . 

In summary, shear flow drives metastable Brownian dis- 



persions to a stationary non-equilibrium state with a multi- 
plicative renormalization of the FDR at long times, which is 
(almost) independent of rescaled time. It nearly agrees for 
variables not advected by flow and takes the universal value 
{jt) = 5 in glasses at small shear rates in leading approxi- 
mation. Corrections arise from the difference of the stationary 
to the transient correlator, and depend on the considered vari- 
able. They alter Xa to values Xa < 5 in the glass. We show 
a new connection between Axa and Ca(i,tiu), see Eq. (|5]), 
which can be tested directly in simulations. 

The derived FDRs characterize the shear-driven relaxation 
at long times, which, according to the ITT approach, is also 
the origin of a (dynamic) yield stress in shear molten glass 
111]. This view captures extended simulations ^ |23tl and 
broad-band experiments 12411 . establishing shear molten glass 
as model for investigating non-equilibrium. Our finding of 
values close to the universal Xqi'yt) = ^ point to intriguing 
connections to critical systems ||9l[l0[]. Open questions con- 
cern establishing such a connection and to address the concept 
of an effective temperature, which was developed for ageing 
and driven mean field models. 

We thank J.-L. Barrat, M. E. Gates and R Ilg for helpful 
discussions. M.K. was supported by the DFG in IRTG 667. 
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This supplementary material summarizes the known ITT 
models for transient correlators used to evaluate the station- 
ary fluctuation functions and susceptibilities presented in the 
main text. 



SCHEMATIC MODEL 

In the schematic -niodel of ITT, a normaUzed transient 
correlator <I>(t) = C^*^ (t) is considered which represents col- 
lective density fluctuations. It obeys the equation of motion 
[12], 



= + / dt'm{j,t-t')^{t'), 

Jq 



t. 



with initialization $(t 0) — > 1 
studied values = 2, = w|(\/4/?;, 
transition at e = 0, and take m(0, t) in order to calculate qui- 
escent (7 = 0) correlators 12511 . For the corresponding sta- 
tionary correlator C{t), we additionally have to estimate the 

Ai) 



We use the much 
1) with the glass 



with Goal f ~ \- Additionally we approximate the normal- 
izations for C^*\t) and C{t), i.e., the respective structure fac- 
tors, to be equal. For incoherent fluctuations, this is exact. 



WAVEVECTOR DEPENDENCE 



To derive the long time FDR as function of wavevector q for 
incoherent, i.e., single particle density fluctuations, we use the 
normalized transient correlator from the ISHSM-model [12]. 
In glassy states, it approximately equals at long times 



CW(t)«/,cxp(^-7^t 



where fq is the g-dependent nonergodicity parameter and hq 
is another coefficient which as well as fq follows from equi- 
librium properties. The quiescent correlator is constant at long 
times, c!f \t) = fq. The long time FDR then follows from 



normalized shear stress a. In flie F^j'^-model, flie shear mod- 7.9^ ^ estimated as above 

ulus is described by the transient correlator and we insert a 
factor of i to account for the different plateau heights of the 
two normalized functions, a is hence estimated as 



/•OO ^ 
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